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C3 ' Abstract 



There is a well developed and useful theory of Hamiltonian reduction for 
semidirect products, which applies to examples such as the heavy top, com- 
pressible fluids and MHD, which are governed by Lie-Poisson type equations. 
In this paper we study the Lagrangian analogue of this process and link it 
with the general theory of Lagrangian reduction; that is the reduction of vari- 
ational principles. These reduced variational principles are interesting in their 
own right since they involve constraints on the allowed variations, analogous to 
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what one finds in the theory of nonholonomic systems with the Lagrange d'Al- 
embert principle. In addition, the abstract theorems about circulation, what 
we call the Kelvin-Noether theorem, are given. 
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1 Introduction 

The main purpose of this paper is to develop the variational structure of the general 
Euler-Poincare equations that are the Lagrangian counterpart to the Lie-Poisson 
equations associated with a semidirect product and to show how this is related to 
the general theory of Lagrangian reduction. We also want to explain the abstract 
notion of circulation that gives rise to a general Kelvin-Noether theorem. 

To accomplish these goals, it will be convenient to first recall some simpler 
situations, namely the two general ways of abstracting the classical Euler equations 
for a fluid or a rigid body, namely Lie-Poisson systems on the dual of a Lie algebra 
and their Lagrangian counterpart, the "pure" Euler-Poincare equations on a Lie 
algebra. 

The Lie-Poisson Equations. The Lie-Poisson equations are Hamiltonian equa- 
tions on the dual of a Lie algebra and represent an abstraction of the Euler equations 
for a rigid body (in body momentum representation) as well as the Euler equations 
for an ideal incompressible fluid (in spatial representation) . This set up and its coun- 
terpart, the Euler-Poincare equations on the Lie algebra are the basic ingredients 
used in the fundamental paper of Arnold [1966a]. 

Ignoring, for simplicity, function space technicalities in the infinite dimensional 
case (see Ebin and Marsden [1970] for one approach to dealing with them), we let 
G be a Lie group with Lie algebra q and let F, K be real valued functions on the 
dual space q* . Denoting elements of q* by fi, let the functional derivative of F at fi 
be the unique element 5F/5fi of q defined by 

lim i[F(/x + e5u) - F(ji)] = (dp, ^ , (1.1) 



for all 5/j, € g*, where (,) denotes the pairing between g* and g. Define the (=L) 
Lie-Poisson brackets by 



{F,G} ± (a)=±(» 



5F <5G 
5/i ' S/j, 



1.2) 



These brackets, discovered by Lie [1890], make g* into a Poisson manifold and 
the Hamiltonian equations associated with a given Hamiltonian H are called the 
Lie-Poisson equations. Thus, the Lie-Poisson equations are determined by F = 
{F, H} for all F. 

We recall these equations in the finite dimensional case in coordinates. Choose 
a basis e±, . . . ,e r of g (so dimg = r). Define, as usual, the structure constants Ci 
of the Lie algebra by [e a , e&] = C^e^, where a, b run from 1 to r and a sum on d 
is understood. If £ € g, its components relative to this basis are denoted £ a , so 
£ = £ a e a . If e 1 , . . . , e n is the corresponding dual basis, and we write [i = fi a e a (with 
a sum understood), the (=L) Lie-Poisson brackets become 

{ F, K}± („) = ± C^|£|I. (1.3) 

The Lie-Poisson equations are 

(l a = {fl a ,H}±(fl) = ±C^ b fl d 



d dH 



da 



b 



or intrinsically, 



A = Tad*8H/6» V- ( L4 ) 



Here ad^ : g — > g is the adjoint map 77 1— > [£, r/], and ad| : g* — > g* is its dual. 

However, the point of view of general Poisson manifolds for these systems is 
sometimes misleading, at least for applications in mechanics. What is more in- 
sightful, and which has its roots in Arnold [1966a] is the point of view that the 
Lie-Poisson brackets arise from canonical brackets on the cotangent bundle T*G by 
reduction. Namely, if we identify g* with the natural Poisson quotient T*G/G using 
the left action of G, we get the minus Lie-Poisson bracket, while the quotient with 
the right action gives the plus Lie-Poisson bracket. 

The rigid body is naturally a left invariant system on T*SO(3), while ideal 
fluids naturally give a right invariant system on the group of volume preserving 
diffeomorphisms of the fluid within the container. 

Poincare and the Euler equations. Poincare's work on the gravitating fluid 
problem continues the line of investigation begun by MacLaurin, Jacobi and Rie- 
mann and was a natural precursor to his famous paper, Poincare [1901], in which he 
laid out the basic equations of Euler type on Lie algebras. He was certainly aware 
that this formalism includes the rigid body, heavy top and ideal fluids as special 
cases. 



To state the Euler-Poincare equations, let q be a given Lie algebra and let I : 
q — > R be a given function (the Lagrangian), let £ be a point in q and let / € Q* be a 
given force (the nature of which we shall explicate later) . Then the Euler-Poincare 
equations are the following equations for the evolution of the basic variable £: 

dSl .,«,. 

The notation is as follows: 51 /5£ £ Q* is the derivative of I with respect to £ — the 
dual to the derivative notation used for the Lie-Poisson equations. 

As with the Lie-Poisson setting, these equations are valid for either finite or 
infinite dimensional Lie algebras. For fluids, Poincare was aware that one needs to 
use infinite dimensional Lie algebras, as is clear from Poincare [1910]. The equations 
also come in two flavors, one for left invariant systems (the one given) and one for 
right invariant systems in which case one changes the sign in front of the adjoint 
operator (or changes conventions in the Lie algebra brackets). 

In the finite dimensional case, the equations read 



d dl _ b dl 
~dtW ~ ad d£ b 



C b ad ^e + fd- (1.5) 



For example, consider the Lie algebra R 3 with the usual vector cross product. (Of 
course, this is the Lie algebra of the proper rotation group in R 3 .) For I : R 3 — > R, 
the Euler-Poincare equations become 

dtdn~ an x " + / ' 

which generalize the Euler equations for rigid body motion to include external forces 
and Lagrangians that needn't be quadratic. 

These equations were written down for a certain class of Lagrangians I by La- 
grange [1788, Volume 2, Equation A on p. 212], while it was Poincare [1901] who 
generalized them (without reference to the ungeometric Lagrange!) to any Lie al- 
gebra. However, Lagrange did grapple (in much of volume 2 of his treatise) with 
the derivation and deeper understanding of the nature of these equations. While 
Poincare may have understood how to derive them from other principles, he did not 
reveal this. 

Of course, there was a lot of mechanics going on in the decades leading up to 
Poincare's time. It is a curious historical fact that the Euler-Poincare equations 
were not pursued extensively and systematically until rather recently. While many 
authors mentioned them (see, e.g., Hamel [1904, 1949] and Chetayev [1941]), it was 
not until the Arnold school that a deeper understanding was achieved and was used 
for purposes of hydrodynamical stability (see Arnold [1966b, 1988] and Arnold and 
Khesin [1992, 1997]). 

We now recall the derivation of the "pure" Euler-Poincare equations (i.e., the 
Euler-Poincare equations with no forcing or advected terms) for left-invariant La- 
grangians on Lie groups (see Marsden and Scheurle [1993a, b], Marsden and Ratiu 
[1994], and Bloch et al. [1996]). The Lagrangian counterpart to the Lie-Poisson 
reduction of Poisson manifolds mentioned above is the following: 



Theorem 1.1 Let G be a Lie group and L : TG — > M a left (right) invariant 
Lagrangian. Let I : q — > R be its restriction to the tangent space at the identity. For 
a curve g(t) E G, let £(£) = g(i) _1 g(i) (respectively £(£) = g{t)g(t)~ l ). Then the 
following are equivalent: 



i Hamilton's principle 



5 / L(g(t),g(t))dt = (1.6) 



<( 



holds, as usual, for variations 5g(t) of g(t) vanishing at the endpoints. 
ii The curve g{t) satisfies the Euler- Lagrange equations for L on G. 

iii The "variational" principle 

b 
l{£{t))dt = (1.7) 

holds on q, using variations of the form 

6S = fi±[Z,rj\, (1.8) 

where n vanishes at the endpoints (+ corresponds to left invariance and — to 
right invariance)!^ 

iv The pure Euler- Poincare equations hold 

d 51 ^ 51 

JtTz - ±ad ^ 



±ad|-. (1.9) 



We make some comments on the proof. First of all, the equivalence of i and 
ii holds on the tangent bundle of any configuration manifold Q, by the general 
Hamilton principle. To see that ii and iv are equivalent, one needs to compute the 
variations 5£, induced on £ = g~ l g = TL g -ig by a variation of g. We will do this 
for matrix groups; see Bloch, Krishnaprasad, Marsden, and Ratiu [1994, 1996] for 
the general case. To calculate this, we need to differentiate g~ l g in the direction of 
a variation 5g. If 5g = dg/de at e = 0, where g is extended to a curve g e , then, 

5£ = — _1 — 
de dt 

while if 7] = g _1 5g, then 



dt de 



d _ 1 d 
dt de 

The difference <5£ — r] is thus the commutator [£,rj\. 

1 Because there are constraints on the variations, this principle is more like a Lagrange d'Al- 
embert principle, which is why we put "variational" in quotes. Of course such problems are not 
literally variational. 



To complete the proof, we show the equivalence of iii and iv. Indeed, using the 
definitions and integrating by parts, 

l(Odt = J —5£dt = J —(f, + ad i ri)dt 

d fSl\ ,jn , 

dt[ji) + ad ^\ vdt 

so the result follows. 

Since the Euler-Lagrange and Hamilton equations on TQ and T*Q are equivalent 
for regular Lagrangians or Hamiltonians, it follows that the Lie-Poisson and Euler- 
Poincare equations are also equivalent under a regularity condition. To see this 
directly, we make the following Legendre transformation from q to Q*: 

Note that 

Sh I 5£\ / 81 6Z\ 

and so it is now clear that the Lie-Poisson and Euler-Poincare equations are equiv- 
alent. However, we shall see shortly that this does not work for examples like the 
heavy top, MHD etc., when we write the equations as Lie-Poisson equations on the 
dual of a semidirect product; the reason is that in these examples, the Hamiltonian 
is degenerate, i.e., if the matrix of second derivatives H Pip . is singular. This is a 
fundamental obstacle and may be viewed as one reason that the Lagrangian version 
of the semidirect product theory is interesting. We shall come to this shortly. 

We close this section by mentioning that many other systems can be put into 
Euler-Poincare form. For example, this can be done for the KdV equations as has 
been shown by Ovsienko and Khesin [1987] (an account may be found in Marsden 
and Ratiu [1994]). The same has been shown for many other shallow water equations 
and equations of continuum mechanics, including the equations of plasma physics 
(cf. Holm, Marsden and Ratiu [1997] and Cendra, Holm, Hoyle and Marsden [1997]). 

Plan of the paper. With the above background in hand, we first set the stage for 
the main goal of the paper, the comparison of the general Euler-Poincare equations 
with Lagrangian reduction, by describing some of the ingredients in more detail. 
The first of these is the Hamiltonian theory of semidirect product reduction in §2. 
The next main ingredient, in §3, is the derivation of the general Euler-Poincare 
equations by the method of reduction of variational principles. In §4 we give the 
Kelvin-Noether theorem; while not essential to the main goal of the paper, it is a 
fundamental result for the Euler-Poincare equations, so it is given for completeness. 
In §5 we describe the general theory of Lagrangian reduction and illustrate the 
methods with Wong's equations and use Lagrangian reduction to give a simple 
proof of the falling cat theorem of Montgomery. Finally, in §6, we show how the 
Euler-Poincare equations are linked with the general theory of Lagrangian reduction. 
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2 Semidirect Product Reduction 

Before getting to the Lagrangian analogue of semidirect product reduction, it is 
useful to summarize the Hamiltonian theory. 

The theory begins with the understanding of various examples, such as the 
heavy top, ideal compressible fluids and MHD (magnetohydrodynamics). Build- 
ing on these examples, the general study of Lie-Poisson equations for systems on 
the dual of a semidirect product Lie algebra was developed by many authors such 
as Sudarshan and Mukunda [1974], Vinogradov and Kupershmidt [1977], Ratiu 
[1980], Guillemin and Sternberg [1980], Ratiu [1981, 1982], Marsden [1982], Mars- 
den, Weinstein, Ratiu, Schmidt and Spencer [1983], Holm and Kupershmidt [1983], 
Kupershmidt and Ratiu [1983], Holmes and Marsden [1983], Marsden, Ratiu and 
Weinstein [1984a,b], Guillemin and Sternberg [1980, 1984], Holm, Marsden, Ratiu 
and Weinstein [1985], Abarbanel, Holm, Marsden, and Ratiu [1986], Leonard and 
Marsden [1997] and Marsden, Misiolek, Perlmutter and Ratiu [1997]. As these and 
related references show, these equations apply to a wide variety of systems such 
as the heavy top, compressible flow, stratified incompressible flow, MHD and the 
dynamics of underwater vehicles. 

Generalities on Semidirect Products. We first recall some definitions. Let V 
be a vector space and assume that the Lie group G acts on the left by linear maps 
on V (and hence G also acts on on the left on its dual space V*). As sets, the 
semidirect product S = G®V is the Cartesian product S = G x V whose group 
multiplication is given by 

{gi,v 1 )(g 2 ,v 2 ) = {gi92,vi +giv 2 ), (2.1) 

where the action of g G G on v G V is denoted simply as gv. The identity element 
is (e,0) where e is the identity in G and inversion is (g, v)^ 1 = (g^ 1 , —g~ 1 v). 

The Lie algebra of S is the semidirect product Lie algebra, s = 0©^, with the 
bracket 

m,vi), (&,!*)] = ([6,6], 6^2 - &vi) , (2.2) 

where we denote actions, such as the induced action of g on V by concatenation, as 
in £it>2- 

The adjoint and the coadjoint actions for semidirect products are given by (see, 
e.g., Marsden, Ratiu and Weinstein [1984a,b]): 

{g,v)(S,u) = (gt,gu-(gS)v). (2.3) 

and 

(g,v){n,a) = (gfi + p* v (ga),ga), (2.4) 

where (g,v) G S = G x V, (£,u) G s = g x V, (fi,a) G s* = g* x V*, g£ = Ad s f, 
gp = Ad* ifj,, ga denotes the induced left action of g on a (the left action of G on 



V induces a left action of G on V * -- the inverse of the transpose of the action on 
V), p v : Q —> V is the linear map given by p v (t;) = £v, and p* : V* — > 0* is its dual. 
For o £ 7*, we shall write, for notational convenience, 

p* v a = v o a G q* , 

which is a bilinear operation in v and a. Continuing to employ the concatenation 
notation for Lie group or algebra actions, the identity 

(rja, v ) = — (v o a , rj) (2-5) 

for all v G V, a G V* and r] G is another way to write the definition of v o a G 0*. 
Using this notation, the coadjoint action reads (5, u)(/x, a) = (<7/x + v o (ga),ga). 

When working with various models of continuum mechanics and plasmas it is 
convenient to work with right representations of G on the vector space V (as in, 
for example, Holm, Marsden and Ratiu [1986]). In this context of course the above 
formalism must be suitably modified. 

Lie-Poisson Brackets and Hamiltonian Vector Fields. For a left represen- 
tation of G on V the ± Lie-Poisson bracket of two functions /, k : s* — > R is given 
by 



{f,k}±(p,a) = ±(/i 



8f_ 8k 
dp,' 5p 



, , 615k 6k 5f 
' ' dp 5a 5p Sa / 



where Sf/Sp G 0, and Sf/5a G V are the functional derivatives of /. The Hamilto- 
nian vector field of h : s* — ► R has the expression 

/ n f 1* 5h Sh \ , . 

X h (p, a) = t I adsh/s^ ~ J^ ° a ' ~ 7^ a ) ( 2 - 7 ) 

Thus, Hamilton's equations (the Lie-Poisson equations) on the dual of a semidirect 
product are given by 

A = Tad^/i±-oa, (2.8) 

d = ± — a, (2.9) 

where overdot denotes time derivative. Again, for right representations of G on V 
the above formulae must be appropriately modified. 

Symplectic Actions by Semidirect Products. We consider a (left) symplectic 
action of S on a symplectic manifold P and assume that this action has an equiv- 
ariant momentum map Jg : P — > 5*. Since V is a (normal) subgroup of S, it also 
acts on P and has a momentum map Jy : P — > V* given by Jy = iy ° J5 , where 
iy : V — > s is the inclusion w h^ (0, u) and iy : s* —> V* is its dual. We think of this 
merely as saying that Jy is the second component of J5. 



We can regard G as a subgroup of 5 by g h- ► (5, 0). Thus, G also has a momentum 
map that is the first component of J5 but this will play a secondary role in what 
follows. On the other hand, equivariance of J 5 under G implies the following relation 
for Jy: 

Jv(gz) = gJv(z) (2.10) 

where we denote the appropriate action of g E G on an element by concatenation, 
as before. To prove ( |2.10| ), one uses the fact that for the coadjoint action of S on 
5* the second component is just the dual of the given action of G on V . 

The Classical Semidirect Product Reduction Theorem. In a number of 
interesting applications such as compressible fluids, the heavy top, MHD, etc., one 
has two symmetry groups that do not commute and commuting reduction by stages 
(a theorem from Marsden and Weinstein [1974]) does not apply. In this more general 
situation, it matters in what order one performs the reduction, which occurs, in 
particular for semidirect products. The main result covering the case of semidirect 
products has a complicated history, some of which has been sketched; we follow the 
version of Marsden, Ratiu and Weinstein [1984a, b] . 

The semidirect product reduction theorem states, roughly speaking, that for the 
semidirect product S = G(§)V where G is a group acting on a vector space V and S 
is the semidirect product, one can first reduce T*S by V and then by G and obtain 
the same result as reducing by S. As above, we let 5 = (S) V denote the Lie algebra 
of S. The precise statement is as follows. 

Theorem 2.1 (Semidirect Product Reduction Theorem.) Let S = G(§)V, 

choose a = (/i, a) E 5* x V* , and reduce T*S by the action of S at a giving the coad- 
joint orbit O a through a E 5* . There is a symplectic diffeomorphism between O a and 
the reduced space obtained by reducing T*G by the subgroup G a (the isotropy of G 
for its action on V* at the point a E V* ) at the point /i a E g* defined by restriction: 
fi> a = MiSaj where Q a is the Lie algebra of G a . 

Reduction by Stages. The preceding result is a special case of a general theorem 
on reduction by stages for semidirect products acting on a symplectic manifold (See 
Marsden, Misiolek, Perlmutter and Ratiu [1997] for this and more general results 
dealing with group extensions and see Leonard and Marsden [1997] for an application 
to underwater vehicle dynamics.) 

As above, consider a symplectic action of S on a symplectic manifold P and 
assume that this action has an equivariant momentum map J5 : P — > 5*. As we 
have explained, the momentum map for the action of V is the map Jy : P —* V* 
given by Jy = i* v o J 5 

We carry out the symplectic reduction of P by S at a regular value a = (fi, a) 
of the momentum map J5 for S in two stages. First, reduce P by V at the value a 
(assume it to be a regular value) to get the reduced space P a := Jy (a)/V . Second, 
form the group G a consisting of elements of G that leave the point a fixed using 
the action of G on V*. One shows (and this step is not trivial) that the group G a 
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acts on P a and has an induced equivariant momentum map J a : P a — ► g*, where g a 
is the Lie algebra of G a , so one can reduce P a at the point fj, a := fi\g a to get the 
reduced space {P a )^ a = ^a l {^a)/{Ga)^ a - 

Theorem 2.2 (Reduction by Stages for Semidirect Products.) The reduced 
space {P a )fi a is symplectically diffeomorphic to the reduced space P a obtained by re- 
ducing P by S at the point a = (/i, a). 

Semidirect Product Reduction of Dynamics. There is a technique for reduc- 
ing dynamics that is associated with the geometry of the semidirect product reduc- 
tion theorem. We start with a Hamiltonian H aQ on T*G that depends parametrically 
on a variable ao £ V* . The Hamiltonian, regarded as a map H : T*G x V* — > M. is 
assumed to be invariant on T*G under the action of G on T*G xV* . This condition 
is equivalent to the invariance of the function H defined on T*S = T*G x V x V* 
extended to be constant in the variable V under the action of the semidirect prod- 
uct. By the semidirect product reduction theorem, the dynamics of H ao reduced by 
G ao , the isotropy group of ao , is symplectically equivalent to Lie-Poisson dynamics 
on s* = q* x V*. This Lie-Poisson dynamics is given by the equations ( [2.8] ) and 
(|2.9| ) for the function h([J,,a) = H{a g ,g^ 1 a) where \i = g^ 1 a g . 

Cotangent bundle reduction. It will be useful to recall a few additional facts 
about cotangent bundle reduction. Following standard notation consistent with 
what we have already used, a symplectic reduced space at momentum value \i is 
denoted P^ := J _1 (/i)/G At . For cotangent bundle reduction we are considering 
P = T*Q and the action of a Lie group G on Q with the standard momentum map. 
The simplest case of cotangent bundle reduction is reduction at the momentum value 
zero in which case one has (T*Q)n=o = T*(Q/G), the latter with the canonical 
symplectic form. If G is abelian, then (T*Q) M = T*{Q/G) but the latter has a 
symplectic structure modified by magnetic terms; that is, by the curvature of the 
mechanical connection. This abelian version of cotangent bundle reduction was 
developed by Smale [1970] and Satzer [1977] and was generalized to the nonabelian 
case in Abraham and Marsden [1978]. Kummer [1981] introduced the interpretations 
of these results in terms of the mechanical connection. This set up was effectively 
used, for example, in Guichardet [1984] and Iwai [1987] to study geometric phases in 
classical molecular dynamics. For additional information on the cotangent bundle 
reduction theorem, see Marsden [1992]. 

When combined with the cotangent bundle reduction theorem, the semidirect 
product reduction theorem is a useful tool. For example this shows directly that the 
generic coadjoint orbits for the Euclidean group are cotangent bundles of spheres 
with the associated coadjoint orbit symplectic structure given by the canonical struc- 
ture plus a magnetic term. In fact, this technique allows one to understand the 
analogue of this geometrical structure for general semidirect products. 

The "bundle picture" begun by these early works was significantly developed 
by Montgomery, Marsden and Ratiu [1984] and by Montgomery [1986], and was 
motivated by work of Weinstein and Sternberg on Wong's equations (the equations 
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for a particle moving in a Yang- Mills field). We shall come to the Lagrangian 
counterpart of this theory shortly. 

3 Lagrangian Semidirect Product Theory 

Introduction. Despite all the activity on the Hamiltonian theory of semidirect 
products, little attention was paid to the corresponding Lagrangian side. Now that 
Lagrangian reduction is maturing (see Marsden and Scheurle [1993a, b]) and useful 
applications are emerging such as to Hamilton's principle asymptotics (see Holm 
[1996]) and to numerical algorithms (Wendlandt and Marsden [1997] and Mars- 
den, Patrick and Shkoller [1997]), it is appropriate to consider the corresponding 
Lagrangian question more deeply. 

The theory is entirely based on variational principles with symmetry. Consistent 
with the variational formulation, note that none of the theorems in this section 
require that the Lagrangian be nondegenerate. The theory is not dependent on the 
previous Hamiltonian formulation, although we shall, of course, make links with it 
under the appropriate regularity assumptions. 

The theorems that follow are modeled after the reduction theorem for the pure 
Euler-Poincare equations given in the introduction, although, as we shall explain, 
they are not literally special cases of it. As in the Hamiltonian case, the main 
distinction between the pure Euler-Poincare equations and the general ones is the 
presence of advected quantities, such as the body representation of the direction of 
gravity in the heavy top and the density in compressible fluids. 

Abstractly, these advected quantities are reflected by the fact that the La- 
grangian L and its reduction /, depend on another parameter a € V* , where, as 
in the Hamiltonian case, V is a representation space for the Lie group G and L has 
an invariance property relative to both arguments. 

As we shall see shortly, the resulting Euler-Poincare equations are not the same 
as the pure Euler-Poincare equations for the semidirect product Lie algebra (§) V* . 

The basic ingredients. We begin with a left representation of Lie group G on 
the vector space V and G acts in the natural way on the left on TG x V*: h(v g , a) = 
(hvg,ha). We assume that we have a L : TG x V* — > R is left G-invariant. In 
particular, if ao € V*, define the Lagrangian L ao : TG — ► R by L aQ (y g ) = L(v g , ao). 
Then L ao is left invariant under the lift to TG of the left action of G ao on G, where 
G ao is the isotropy group of ao- 

Left G-invariance of L permits us to define the reduced Lagrangian I : q x V* — ► M 
by l(g~ 1 v g , g~ l a) = L(v g , a). Conversely, this relation defines for any I : q x V* — ► M. 
a left G-invariant function L : TG x V* —> R. 

For a curve g(t) € G, let £(£) := g{t)~ l g{t) and define the curve a(t) as the 
unique solution of the following linear differential equation with time dependent 
coefficients a(t) = —£(t)a(t), with initial condition a(0) = ao- The solution can be 
written as a(t) = <7(t) _1 ao- 
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The Euler-Poincare equations. The generalization of the pure Euler-Poincare 
theorem is the following. 

Theorem 3.1 With the preceding notation, the following are equivalent: 

i With ao held fixed, Hamilton's variational principle 

5 I" L ao (g(t),g(t))dt = (3.1) 

holds, for variations 5g(t) of g(t) vanishing at the endpoints. 
ii g(t) satisfies the Euler-Lagrange equations for L aQ on G. 
iii The constrained variational principle (of Lagrange d'Alembert type) 

5 I* l(£(t),a(t))dt = (3.2) 

Jtx 

holds on q x V* , using variations of £ and a of the form 

5£ = V + [€, V], $a = ~W, ( 3 - 3 ) 

where n(t) G Q vanishes at the endpoints. 

iv The Euler Poincare equations hold on g x V* 

d 51 ,*5l 51 , ... ,.,. . . .„ ,. 

— — = ad|— + — oo, where a(t) = -£(t)a(t). (3.4) 

at o£ ^ oj; da 

Proof. The equivalence of i and ii holds for any configuration manifold and so, 
in particular, it holds in this case. 

Next we show the equivalence of iii and iv. Indeed, using the definitions, inte- 
grating by parts, and taking into account that rj(t\) = 77(^2) = 0, we compute the 
variation of the integral to be 

S J\{i{t),a(t))dt = f"(/^,S(\ + /sa,^\'\dt 

= r((i'" +ads ")"(" a '^)) <ft 
- r«^(i) + ^i'"> + (i»°'"» dt 

f t2 I d ( 5l\ „<H 51 \ , 



and so the result follows. 

Finally we show that i and iii are equivalent. First note that the G-invariance 
of L : TG xP-tR and the definition of a{t) = g(t)~ 1 ao imply that the integrands 
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in (3.1) and ([T^) are equal. However, all variations 5g(t) £ TG of g(t) with fixed 



endpoints induce and are induced by variations 5£(t) € (j of £(£) of the form 5£ 
V + [£> v] with r/(i) £ g vanishing at the endpoints; the relation between Sg(t) and 
r)(t) is given by r)(t) = g(t)~ 1 5g(t). This is the content of the following lemma 
(which is elementary for matrix groups) proved in Bloch et al. [1996]. 

Lemma 3.2 Let g : U Cl 2 — > G be a smooth map and denote its partial derivatives 
by £(jk, e) = TL 9(t>£) -i (dg{t, e)/dt) and n(t, e) = TL g ^-i (dg(t, e)/de). Then 

Conversely, if U is simply connected and £, rj : U — > q are smooth functions sat- 
isfying HjS.dj) then there exists a smooth function g : U — > G such that £(£, e) = 
TL g{t£) -i(dg(t,e)/dt) andn(t,e) = TL g{t£) -i(dg(t,e)/de). 

Thus, if i holds, we define n(t) = g(t)~ 1 5g(t) for a variation Sg(t) with fixed 
endpoints. Then if we let <5£ = fi'(i) _1 5(^)) we have by the above proposition 5£ = f]+ 
[£,rj\. In addition, the variation of a(t) = g(t) oq is Sa(t) = —n(t)a(t). Conversely, 
if 5£ = rj + [£, n] with 77(f) vanishing at the endpoints, we define Sg(t) = g(t)rj(t) 
and the above proposition guarantees then that this 5g(t) is the general variation of 
g(t) vanishing at the endpoints. From Sa(t) = —r](t)a(t) it follows that the variation 
of g(t)a(t) = a® vanishes, which is consistent with the dependence of L ao only on 

g(t),g(t). ■ 

Remark. The Euler-Poincare equations are not the pure Euler-Poincare equations 
because we are not regarding q x V* as a Lie algebra. Rather these equations 
are thought of as a generalization of the classical Euler-Poisson equations for a 
heavy top, but written in body angular velocity variables. Some authors may thus 
prefer the term Euler-Poisson- Poincare equations for these equations. The following 
argument shows that these Euler-Poincare equations (|3.4| ) are not the pure Euler- 
Poincare equations for the semidirect product Lie algebra 0©V*. Indeed, by ( |1.9[ ) 
the pure Euler-Poincare equations 

d & 1* M fj. \ ^Tr* 

dtW^) =ad ^W^)> (e < a)GS ® y 

for I : g (J) V* — ► R become 

d 51 „5l 81 d 51 51 

= ad f 1 o a, = — £ — , 

dt o£ ? <5£ 5a dt 5a 5a 

which is a different system from that given by the Euler-Poincare equation ( |3.4[) 
and a = — £a, even though the first equations of both systems are identical. 
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The Legendre Transformation. As we explained earlier, one normally thinks of 
passing from Euler-Poincare equations on a Lie algebra g to Lie-Poisson equations 
on the dual q* by means of the Legendre transformation. In our case, we start 
with a Lagrangian on q x V* and perform a partial Legendre transformation in the 
variable £ only, by writing 

^ = 77' /l (/ i ' a ) = (m>0 -J(£>°0- ( 3 - 6 ) 

A simple calculation shows that 5h/5fi = £, and 5h/5a = —51/ 5a, so that ( [3.4] ) and 
d(i) = — £(i)a(i) imply (2.7) for the minus Lie-Poisson bracket (that is, the sign 



+ in ( |2.7| )). If this Legendre transformation is invertible, then we can also pass 
from the minus Lie-Poisson equations ( |2.7D to the Euler-Poincare equations ( |3.4[ ) 
together with the equations a{t) = —£(t)a(t). 

4 The Kelvin-Noether Theorem 

In this section, we derive a version of the Noether theorem that holds for solutions 
of the Euler-Poincare equations. Our formulation is motivated and designed for 
continuum theories (and hence the name Kelvin-Noether), but it holds in a more 
general situation. Of course it is well known (going back at least to the work of 
Arnold [1966a]) that the Kelvin circulation theorem for ideal flow is closely related 
to the Noether theorem applied to continua using the particle relabeling symmetry 
group. 

The Kelvin-Noether Quantity. Start with a Lagrangian L ao depending on a 
parameter oq € V* as above. Consider a manifold C on which G acts (we assume 
this is also a left action) and suppose we have an equivariant map /C : C x V* — ► q**. 

In the case of continuum theories, the space C may be chosen to be a loop space 
and (IC(c, a) , fj) for c E C and /i€g* will be a circulation. This class of examples 
also shows why we do not want to identify the double dual q** with q. 

Define the Kelvin-Noether quantity I : C x g x V* — >IRby 

I(c,H,a) = (}C(c,a),^,a)y (4.1) 

Theorem 4.1 (Kelvin-Noether) Fixing cq E C, let £(t),a(t) satisfy the Euler- 
Poincare equations and define g(t) to be the solution of g(t) = g(t)£(t) and, say, 
g(0) = e. Let c(t) = ^(t)" 1 ^ and I(t) = I(c(t),£(t),a(t)). Then 

j t I(t) = (K(c(t),a(t)),^oay (4.2) 

Proof. First of all, write a(t) = g(t)~ 1 ao as we did previously and use equivariance 
to write I(t) as follows: 



K{c{t), a(t)), ^(£(i), a(t)) ) = ( £(c , a ), g(t) 
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The g^ 1 pulls over to the right side as g (and not g~ l ) because of our conventions 
of always using left representations. We now differentiate the right hand side of this 
equation. To do so, we use the following well known formula for differentiating the 
coadjoint action (see Marsden and Ratiu [1994], page 276): 



d_ 
df 



\g(t)n(t)] = g(t) 



ad £mM*) + -jiM*) 



*«*) 



dt 



where, as usual, 

m=9(t)- 1 g(t). 
Using this and the Euler-Poincare equations, we get 



dt 



— (JC(co,a ),g(t) 



fc{co,a ), — \g(t) 



|«(*),a(*)) 



,* SI „ 51 51 



5l_ 
6a 



o a 



JC(c ,a ),g(t) 
fZ(co,a ),g(t) 

5(t) _1 /C(c ,a 
/C(c(t),a(t)) 



where, in the last steps we used the definitions of the coadjoint action, the Euler- 
Poincare equation ( |3.4| ) and the equivariance of the map /C. ■ 





'51 




)), 


— o a 


5a 


'51 


\ 


— Offl 


) 


5a 
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Corollary 4.2 For the pure Euler-Poincare equations, the Kelvin quantity I(t), 
defined the same way as above but with I : C x g — ► R, is conserved. 

The Heavy Top. As a simple illustration of the results so far, we consider the 
heavy top. For continuum examples, we refer to Holm, Marsden and Ratiu [1997] for 
fluid theories and to Cendra, Holm, Hoyle and Marsden [1997] for Vlasov plasmas. 
The heavy top kinetic energy is given by the left invariant metric on SO(3) whose 
value at the identity is (x, y) = Ix • y, where x, y € R 3 are thought of as elements 
of so (3), the Lie algebra of 50(3), via the isomorphism x G R 3 i— > x € so (3), 
xy := x x y, and where I is the (time independent) moment of inertia tensor in body 
coordinates, usually taken as a diagonal matrix by choosing the body coordinate 
system to be a principal axes body frame. This kinetic energy is thus left invariant 
under the full group 50(3). The potential energy is given by the work done in lifting 
the weight of the body to the height of its center of mass, with the direction of gravity 
pointing downwards. If m denotes the total mass of the top, g the magnitude of the 
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gravitational acceleration, % the unit vector of the oriented line segment pointing 
from the fixed point about which the top rotates (the origin of a spatial coordinate 
system) to the center of mass of the body, and I its length, then the potential energy 
is given by — mglH~ 1 e3 ■ %, where e3 is the axis of the spatial coordinate system 
parallel to the direction of gravity but pointing upwards and where R E SO(3) is 
the orientation of the body. This potential energy breaks the full SO (3) symmetry 
and is invariant only under the rotations S 1 about the e3~axis. 



To apply Theorem 3.1, we take the Lagrangian of the heavy top to be the 
kinetic minus the potential energy, regarded as a function on TSO(3) x R 3 -> K. In 
particular, the potential is given by £/(i4R, v) = mgl'R^ 1 v • %, where ur. G TrSO(3) 
(eventually identified with R) and v £ R 3 (eventually identified with the direction of 
gravity), which is easily seen to be S r O(3)-invariant. Thus, the heavy top equations 
of motion in the body representation are given by the Euler-Poincare equations 
(|3.4| ) for the Lagrangian I : so (3) x R 3 — > R. To compute the explicit expression of 
I, denote by ft the angular velocity and by II = Ifi the angular momentum in the 
body representation. Let T = R _1 v; if v = e3, the unit vector pointing upwards on 
the vertical spatial axis, then T is this unit vector viewed by an observer fixed and 
moving with the body. The Lagrangian I : so (3) x R 3 — ► R is thus given by 

J(0, r) = L(R" 1 n R , R _1 v) = -n fl + mglT ■ x ■ 

It is now straightforward to compute the Euler-Poincare equations. First note that 

SI „ SI 
- = n, - = mgl X . 

Since 

ad^ n = nxfi, vor = -Txv,and QT = -T x Q , 

the Euler-Poincare equations are 

II = II x fj + mglT x x 

which are coupled to the T evolution T = T x fi . This system of two vector equations 
for IT, r are the classical Euler-Poisson equations and they describe the motion of 
the heavy top in the body representation. 

To illustrate the Kelvin-Noether theorem, choose C = q and let /C : C x V* — > 
q** = q be the map (W, T) h-> W. Then the Kelvin-Noether theorem gives the 
statement 

j t (W,U) = mgl(W,Tx X ) 

where W(t) = R(t) _1 '\v; in other words, W(t) is the body representation of a 
space fixed vector. This statement is easily verified directly. Also, note that 
(W, II) = (w,7r), with tv = R(i)II, so the Kelvin-Noether theorem may be viewed 
as a statement about the rate of change of the momentum map of the heavy top 
(that is, the spatial angular momentum) relative to the full group of rotations, not 
just those about the vertical axis. 
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5 General Lagrangian Reduction. 

The Lagrangian analogue of cotangent bundle reduction was developed by Cendra, 
Ibort and Marsden [1987] and Marsden and Scheurle [1993a, b]. (Routh, already 
around 1860 investigated in coordinates what we would call today the abelian ver- 
sion.) These references developed the coordinate version of the resulting reduced 
Euler-Lagrange equations as well as the corresponding associated reduction of vari- 
ational principles. Another interesting paper in this subject is that of Weinstein 
[1996] who used groupoids to show, amongst other things, one can regard discrete 
and continuous Lagrangian reduction from a single viewpoint. Discrete Lagrangian 
systems with symmetry, following ideas of Veselov and Moser are of great interest in 
numerical algorithms (see Wendlandt and Marsden [1997] and references therein), 

One starts with a G- invariant Lagrangian L on TQ, which induces a Lagrangian 
I on the quotient space (TQ)/G, An important point, but which is easy to see, is 
that, assuming the group actions are free and proper, this quotient space (TQ)/G 
is intrinsically a bundle over T(Q/G) with a fiber modeled on the Lie algebra g. 
Denote, in a local trivialization, the variables in the base, or shape space Q/G by 
r a 7 r a and the fiber variables by £ a . In such a local trivialization, the reduced 
equations are the Hamel equations: 

d dl dl 



dt dr a dr a 

±^_ d J_ C a pd n 

dtoe de 

However, for applications to stability problems as well as being global, it is 
desirable to choose a connection A, say the mechanical connection, on the bundle 
Q ~* Q/G. Then the variables in the shape space and the vertical part of a 
velocity vector on Q become globally and intrinsically defined. In coordinates, the 
vertical part is given by Vl a = £ a + A^(r)f a where the components of the connection 
with respect to the local coordinates are denoted A a a . In such variables, the resulting 
reduced Euler-Lagrange equations are 



d dl dl dl 

dtdjz ~ dr^ 'dl^ K ~ 1Jaf3 



- ( R a r? + F a Q d \ 



where B^p is the curvature of the connection A a , C^ d are the structure constants 
of the Lie algebra g and where £® d = C^ d A b a . This second set of equations are 
generalizations of the Euler-Poincare equations. The two sets of equations are called, 
respectively, the horizontal and vertical equations. 

This theory has had a large impact on, for example, the theory of nonholonomic 
systems (see Bloch, Krishnaprasad, Marsden and Murray [1996]). It has also proven 
very useful in optimal control problems. For example, it was used in Koon and 
Marsden [1997a] to extend the falling cat theorem of Montgomery [1990] to the case 
of nonholonomic systems. 
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The work of Cendra, Marsden and Ratiu [1997] develops the geometry and varia- 
tional structure of the reduced Euler-Lagrange equations. This work starts with the 
following view of the bundle picture. As above, choose a connection A on Q — > Q/G 
such as the mechanical connection and let q denote the vector bundle over shape 
space Q/G which is the associated bundle to q and the adjoint action of G on q. 
There is a bundle isomorphism 

TQ/G^T(Q/G)@~q 

determined by the projection on the first factor and by the vertical projection asso- 
ciated with the connection on the second. The sum is a Whitney sum (i.e., fiberwise 
a direct sum) of vector bundles over Q/G. 

Using the geometry of this bundle TQ/G = T(Q/G) © 0, one obtains a nice 
intrinsic interpretation of the above reduced Euler-Lagrange equations in terms of 
covariant derivatives of induced connections. One easily gets the dynamics of par- 
ticles in a Yang-Mills fields as well as many other interesting examples (such as a 
rigid body with rotors, etc.) as special cases. We will illustrate this with Wong's 
equations below. Preliminary calculations show that in nonholonomic mechanics we 
will get a beautiful geometric interpretation of the important momentum equation 
of Bloch, Krishnaprasad, Marsden and Murray [1996]. 

The dual bundle picture. The above construction sheds light on the bundle 
picture for cotangent bundles mentioned earlier. Taking the dual of the above 
isomorphism gives an isomorphism 

T*Q/G^T*(Q/G)®g*. 

Like any quotient Poisson manifold, the space (T*Q)/G has a natural Poisson struc- 
ture. The description of this Poisson structure viewed on the bundle T*(Q/G) © 5* 
involves a synthesis of the canonical bracket, the Lie- Poisson bracket, and curvature 
and is a very interesting way of encoding the original description of Montgomery, 
Marsden and Ratiu [1984] of the Poisson structure on (T*Q)/G. One gets from 
this picture both the usual Lie- Poisson description of dynamics and the Hamiltonian 
description of semidirect product theory. 

Reduced variational principles. The above picture also leads naturally to a 
generalization of the reduced variational principles described for the Euler-Poincare 
equations above. One does this by dividing the variations in Hamilton's principle 
on T*Q into horizontal and vertical parts, which then drop appropriately to the 
quotient space in a manner similar to that in the Euler-Poincare case. 

We regard this reduction of variational principles as fundamental. Another ap- 
proach to reduction of Lagrangian systems is based on reducing geometric objects 
such as almost tangent structures (see for example, De Leon, Mello and Rodrigues 
[1992] and subsequent papers of this group). The present framework seems to fur- 
ther clarify this, and identifies exactly which part of the reduced space is a tangent 
bundle and which part is of Euler-Poincare type. 
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Stability under reduction. One of the motivations is to produce a category, a 
Lagrangian analogue of the symplectic and Poisson category, that is stable under 
reduction. As the Euler-Poincare equations show, this category cannot be that of 
tangent bundles or even locally tangent bundledike objects. 

However, the above ideas lead to an answer to this question in a natural way. We 
enlarge the traditional category for Lagrangian mechanics, namely tangent bundles, 
to bundles of the form 

TQ(BU, 

where U is a vector bundle with a connection over Q, each fiber of which carries a Lie 
algebra structure and where the base Q also carries a two form (a "magnetic term"). 
This structure is rich enough to include the reduced Euler-Lagrange equations (by 
taking U = g) and carry variational methods. 

In addition, this category is stable under reduction. If one reduces this structure 
by the action of a group K, then Q gets replaced by Q/K, while the fiber expands 
by enlarging each Lie algebra fiber by a group extension construction, with the base 
two form playing the role of the cocycle. This provides a Lagrangian analogue of the 
reduction by stages procedure of Marsden, Misiolek, Perlmutter and Ratiu [1997]. 

Wong's Equations. As another illustration of the ideas of Lagrangian reduction, 
we consider the important example of Wong's equations. See Montgomery [1984] 
and the references therein for the original papers. This example is important for its 
own sake, but also through the fact that it enters into some fundamental optimal 
control problems, as shown by Montgomery [1990] and Koon and Marsden [1997a]. 

We begin with a Riemannian manifold Q with a free and proper isometric action 
of a Lie group G on Q. Let A denote the mechanical connection; that is, it is the 
principal connection whose horizontal space is the metric orthogonal to the group 
orbits. The quotient space Q/G = X inherits a Riemannian metric from that on Q. 
Given a curve c{t) in Q, we shall denote the corresponding curve in the base space 
X by r(t). 

The optimal control problem under consideration is as follows: 

Isoholonomic Problem (Falling Cat Problem). Fixing two points q±, #2 G Q, 
among all curves q(t) € Q, < t < 1 such that q(0) = qo, <?(1) = qi and q(t) € hor g ( t ) 
(horizontal with respect to the mechanical connection A), find the curve or curves 
q(t) such that the energy of the base space curve, namely, 

1 rl 



, \\ r \\ dt, 

1 ./0 



is minimized. 



Theorem 5.1 (Montgomery (1990, 1991).) If q(t) is a (regular) optimal tra- 
jectory for the isoholonomic problem, then there exists a curve X(t) E g* such that 
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the reduced curve r(t) in X = Q/G together with X(t) satisfies Wong's equations: 



Pa 

A 6 



A B a r ldaP \,m 



\ (~ia Ad -a 

where g a p is the local representation of the metric on the base space X; that is 

2 



1 



1 



\r\ 



gapr a r p , 



2" " T 
g 01 is the inverse of the matrix g a p, p a is defined by 

•J 



01 



Pa = 0r^ = 9afi 
and where we write the components of A as A b a and similarly for its curvature B. 

Proof. By general principles in the calculus of variations, given an optimal solution 
q(t), there is a Lagrange multiplier \(t) such that the new action function defined 
on the space of curves with fixed endpoints by 



©[?(■)] 



l\\r(t)\\ 2 + (\(t),Aq(t)) 



dt 



has a critical point at this curve. Using the integrand as a Lagrangian, identify- 
ing f2 = Aq and applying the reduced Euler-Lagrange equations to the reduced 
Lagrangian 

l(r,r,tt) = ^\\r\\ 2 + (\,n) 
then gives Wong's equations by the following simple calculations: 



dl 
0r a 



g a pr 



P. 



01 



1 dgfr . 



rY*t-* <y* y • 



01 



A„. 



0r a 2 0r a ' QQ, a 

The constraints are Q = and so the reduced Euler- Lagrange equations become 
d 01 01 



dt 0r a 0r a 



dt 



-MbV) 



\ (~ta Ad -a 
A a { ^db JK a r ■ 



But 



d 01 



01 



dt 0r a 0r a 



ldgfa 

P° 2 0r« r r 



Pa + 77 



2 0r c 



2 0r c 



-PpP~n 
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and so we have the desired equations. ■ 

Let us now give another version of this procedure in abstract notation and, at the 
same time, show how it fits in with the bundle view of Lagrangian reduction. Again, 
let it : Q — > X be a principal bundle with structure group G, a Lie group acting 
on the left, let A be a principal connection on Q and let B be the curvature of A. 
Suppose that g is a given Riemannian metric on X and let V be the corresponding 
Levi-Civita connection. Assume, for simplicity, that G is a compact group with a 
bi-invariant Riemannian metric K. 

Define the Lagrangian L : TQ — > R by 

L(q, q) = X -K (A(q, q), A(q, q)) + X -g (*(?)) (Tir(q, q),Tir(q, q)) . 

An element of 5 has the form v = [q, v]q where q £ Q an d w£8 and where [<?, v]g 
denotes the equivalence class of the pair (g, v) with respect to the group action by 
G. Since K is bi-invariant, its restriction to q is Ad-invariant, and therefore we can 
define K ([q,v] G , [q, v] G ) = K(v,v). 
The reduced bundle is 

T(Q/G) © g = TX g 

and a typical element of it is denoted (x,x,v). The reduced Lagrangian is given by 

l(x,x,v) = -K(v,v) + -g(x)(x,x) . 

Now we will write the vertical and horizontal reduced Euler-Lagrange equations. 
An arbitrary variation 8v in the direction of the fiber of q is of type 8v = [q, 8v]q, 
where Sv 6 is arbitrary. We have 

— (x, x, v)5v = K(v, 5v) 
ov 

and 

dl _ \ 

adi— -(x, x, v) ) 5v = K (v, [v, Sv]) . 
ov J 

Since K is bi-invariant, we have K(ad w u, v) + K(u, ad w v) = 0, for all u, v, and 
w in g. Therefore, K (v,[v,5v]) = —K([v,v],Sv) = 0. Thus the vertical reduced 
Euler-Lagrange equation is one of Wong's equations 

!*<„,. > = a 

This actually agrees with the previous form of Wong's equations since A represents 
the value of the body momentum while the form K(v, • ) represents the spatial mo- 
mentum, which is conserved. After some straightforward calculations, the horizontal 
reduced Euler-Lagrange equation becomes in this case 



(V^) b = ^(^0r)(i, •)) 



which is the other of Wong's equations. Here B denotes the curvature thought of as 
a Lie algebra valued two form on the base. 
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6 The Euler-Poincare Equations via Lagrangian Reduc- 
tion 

Now we are ready to explain how the Euler-Poincare equations can be derived using 
the techniques of Lagrangian reduction. This should be plausible from what we 
have already said simply because the coordinate expression for the reduced Euler- 
Lagrange equations include a generalized form of Euler-Poincare equations. But we 
need to make this connection precise. 

The set up. We consider a slight generalization of the Euler-Poincare equations 
given earlier. This time we consider a group G acting on a vector space V and hence 
on its dual V*. We also consider a configuration space Q, but in this section, G will 
act trivially on Q. Consider a Lagrangian 

L : T{G x Q) x V* -> R 

where G is a group, Q is a manifold and V* is the dual of the vector space V. The 
value of L at the point (g, q, g, q, ao) £ T{G xQ)xV* will be denoted L(g, q, g, q, ao), 
as usual, and, as in the discussion of the Euler-Poincare equations earlier, we think 
of ao as being a parameter that remains fixed along the evolution of the system. 
By construction, the action of G on V and V* satisfies the property 

(ga,gv) = (a,v) 

for all a £ V* , all v £ V and all g £ G. (We will often write (a,v) = (v,a)). 
Assume that L has the following invariance property: 

L(g'g, q, g'g, q, g'a ) = L(g, q, g, q, o ) 

for all ao £ V*, all q £ Q and all </, g £ G. Let 

L(e,q,£,q,a) = l(£,q,q,a), 

for all £ £ q, all q £ Q and all a £ V*. Then the invariance property implies, for all 
g £ G, all q £ Q and all a £ V*, 

L(g,q,g,<i,ao) =l(t,q,q,a) 

where £ = g~ l g and a = g _1 ao. 

A small generalization of the argument given earlier for the Euler-Poincare equa- 
tions proves the next result. 

Theorem 6.1 The following conditions are equivalent: 

(i) The curve (g(t),q(t),ao) is a critical point of the action 

L(g,q,g,q,ao)dt 

to 

with restrictions on variations given by 5g(ti) = for i = 0,1, Sq(ti) = for 
i = 0,1, and Sao = 0. 
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(ii) The curve (£(t) , q(t) , a(t)) , where a(t) = g 1 (t)ao for allt and £(t) = g 1 (t)g, 
is a critical point of the action 

K€,Q,Q,a)dt 

to 

with restrictions on variations given by 

5£ = rj+ [£, rj\ 
where n is any curve on q such that n(ti) = for i = 0, 1, 

Sq(U) = 
for i = 0, 1 and 

5a = —na 

and, besides, the curve a(t) must satisfy 

a + £a = 

for all t. (As earlier, this last condition comes from the condition do = 
together with ao = ga.) 

A direct application of (ii) leads to the reduced equations 

d 51 ,„ SI 51 

together with the standard Euler-Lagrange equations for q: 

dl d dl _ 
dq dt dq 

Recall that, by definition, 

(voa,C) = -((a,v) 

for all £ £ 0, all a £ y* and all w £ V. The first of these equations is of course the 
Euler-Poincare equation. These and the Euler-Lagrange equations for q, together 
with the equation a + £,a = form a complete set of equations of the system in terms 
of the variables (£,<?, a). This framework is applied to Vlasov plasmas in Cendra, 
Holm, Hoyle and Marsden [1997]. 

Now we shall recast conditions (i) and (ii) into an equivalent form. The idea is to 
introduce the condition that ao is conserved by making it the momentum conjugate 
to a cyclic variable, just as one does for the charge in Kaluza-Klein theory. Thus, 
let us define the Lagrangian 

L : T(G x Q x V* x V) -> R 
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by 

L(g, Q, ao,v ,9, q, a , v ) = L{g, q, g, q, a ) + (a , v ) . 

Notice that for the Lagrangian L, ao is the momentum conjugate to the cyclic 
variable Vq, and hence is a constant. The variable Vq is not constant, but its evolution 
is the first order equation vq + dL/dao = 0. Thus, the Euler-Lagrange equations for 
L are equivalent to the Euler-Lagrange equations for L with the parameter ao fixed 
(together with the equation for vo). 

Thus, the strategy is to perform reduction of L by the equivalent process of 
performing standard tangent bundle Lagrangian reduction of L. 

Now we observe that 

GxQxV* xV ^{GxQxV* x V)/G 

is a principal bundle with structure group G acting as before, that is, g'(g, q, oq,vo) = 
(g'g,q,g'ao,g'vo). Moreover, this bundle is isomorphic, as a principal bundle, to the 
trivial bundle 

GxQxV*xV^QxV*xV 

where the action of G is given by g' ■ (g,q,a,v) = (g'g,q,a,v), for all g',g G G, all 
a G V* and all v £ V. This assertion follows using the isomorphism tp : G x Q x 
V* xV -> GxQxV* xV given by t/j(g,q,a ,v ) = (ssg^^acff" 1 ^) = (g,q,a,v). 
It is easy to check that tp (g'(g,q,a ,v )) = g' ■ ip(g,q,a ,v ) = g' ■ (g,q,a,v), for 
all g',g € G, all ao £ V* and all vq G V. One also checks that the composition 
L o Ttp^ 1 =: L is given by 

L V (g, q, a, v, g, q, d, v) = L(g, q, g, q, a) + (a,v + g' l gv) . 

/,From now on we will use the trivial bundle described above and the Lagrangian 

L v . 

Remark. Using techniques like those in Cendra and Marsden [1987], which give 

a version of the Lagrange multiplier theorem, we can show that conditions (i) and 

(ii) are equivalent to either of the following conditions, 

(in) The curve (g(t),q(t),ao,vo) is a critical point of the action 

rti _ 

L(g, q,a , v , g,q, d , v )dt 
to 

for variations satisfying the endpoint conditions 5g(ti) = , Sq(ti) = , 
5ao(ti) = and Svo(U) = for i = 0,1. 

(iv) The curve (g(t),q(t),a(t),v(t)) is a critical point of the action 

L v (g,q,a,v,g,q,a,v)dt 

to 

for variations satisfying the endpoint conditions Sg(ti) = 0, Sq(ti) = 0, 8a(U) = 
and5v(ti) = fori = 0,1. ♦ 
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We shall show that the reduced Euler-Lagrange equations for the Lagrangian 
L v are the same as the reduced system for I. To do this, we first calculate the 
reduced bundle T(Q x V* x V) © q. Since the bundle G x Q x V* x V is trivial, 
the bundle q is q = Q x V* x V x g. The Lie algebra structure on g is given by 

[(q,a,v,£i),(q,a,v,£ 2 )] = (q,a,,v,[ti,&])- 

Let us choose the trivial principal connection A on G x Q x V* x V, that is, the 
connection given by A(g, q, a, v, g, q, a, v) = g~ x 'g. Since the connection is trivial, it 
has zero curvature. The reduced Lagrangian is clearly given by 

l v (£, q, a, v, q, a, v) = /(£, q, q, a) + (a,v + £v) . 

Now one is ready to calculate the reduced Euler-Lagrange equations for l v using 
this set up. This can be done in coordinates using the coordinate version of the 
reduced Euler-Lagrange equations we gave earlier, or in their intrinsic formulation. 

Since the connection is trivial and its curvature is zero, the horizontal equations 
for / are simply the usual Euler-Lagrange equations for these variables; these give 
the equations for q as well as the equations for a and v. The Euler-Lagrange equation 
for q with respect to / is clearly the same as the Euler-Lagrange equation for q 
with respect to I. The horizontal equation for a is the equation 

d dl v dl v 



dt da da 

which, since I does not depend on d, is simply 

■ , dl 
v + £v + — = 0, 
oa 

which is equivalent to the equation ^o + dL/dao = that we had before. 
The Euler-Lagrange equation for v is 

d dl v dl v 



dt dv dv 

Since dl v /dv = a and dl v /dv = — £a, this gives the correct equation a + £a = for 
a. 

The vertical equation becomes 

d dl v _ (dl v 

dt~bj ~ ad « V^ 

Clearly 

dl v dl 

-W = di + voa > 



and so the vertical equation becomes 



d dl . . dl 1+ . . 

— —-+voa + voa = adt— - + achmoa) 
dt d£ ? dt, 
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Using the equations for d and v, this becomes 

i.e., 

d dl ^dl dl . 

di~dl ~ Z~dl ~ d~ ° ° = ^ ' oa + v<> ^ a > + ad ^ v ° a >- 

But the right hand side is identically zero, as is seen from the definitions. Thus, we 
have achieved our goal of recovering the Euler-Poincare equations with an advected 
parameter from Lagrangian reduction. 
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